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Abstract 


A  three-level  system  of  intermediate  transition  frequencies 
“^2  end  “23*  with  nonvanishing  dipole  matrix  elements  for  these 
transitions  only,  is  considered  to  be  coherently  pumped  at  a 
frequency  near  1/2(U12  *  “23>  end  radiating  at  both  intermediate 
frequencies.  A  sharply  resonant  rise  in  the  second-  and 
third-level  population  is  found  for  slight  detuning  of  the  pump 
from  1/2(U12  *  “23>  •  Population  inversion  in  either  transition 
is  shown  to  be  possible,  depending  on  the  pertinent  parameters. 
The  induced  dipole  moment  as  well  as  the  coherence  of  the 
radiation  is  discussed. 
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INTRODUCTION 


Production  of  population  inversion  by  means  of  a 
“three-level**  scheme  is  one  of  the  standard  methods  used  in 
lasers  and  masers  pumped  by  an  electromagnetic  field.  In  this 
scheme,  the  <1,3>  transition  and  either  the  <1,2>  ox.  the  (2,3) 
transition  couple  to  the  electromagnetic  field,  while  the  other 
intermediate  transition  is  nonradiative  and  must  be  sub3ected  to 
a  perturbation  that  produces  relaxation.  The  pump  acts  on  the 
<1.3)  transition,  and  a  steady-state  population  at  the 
intermediate  radiative  frequency  results  for  an  appropriate  range 
of  the  pertinent  parameters,  as  follows  from  conventional  rate 
equations.  In  order  that  this  scheme  be  applicable,  the  states 
corresponding  to  the  first  and  third  level  must  be  of  opposite 
parity,  and  the  pump  frequency  must  be  higher  than  the  laser 
frequency.  In  the  present  report  an  entirely  different  method  of 
steady-state  population  inversion  in  a  three-level  system  will  be 
described,  applicable  when  the  top  and  bottom  states  are  of  the 
same  parity,  both  intermediate  transitions  are  radiative,  and  the 
pump  frequency  may  be  either  higher  or  lower  than  the  frequency 
for  which  inversion  is  produced. 

II.  FORMULATION  OF  PROBLEM  AND  DERIVATION  OF  EQUATIONS  OF  MOTION 
We  consider  an  atomic  system  ("atom")  of  three  energy 
levels,  'ftui,  -hu)2*  -H  w3  in  ascending  order,  with  relatively  close 
intermediate  transition  frequencies  ^3*  <wi3  -  t  I  ) 

and  only  two  nonvanishing  dipole  matrix  elements  ond  V23* 

corresponding  to  the  two  intermediate  transitions,  respectively. 

A  pump  field  E  *  2Eq  cosut,  with  I  w  -  <  1  /2)  u  1 3  I  <  <  <*>,  acts  on  the 


atom,  which.  In  turn,  ia  either  coupled  to  two  loaay  cavity  modes 
with  resonant  frequencies  wi2  and  “>23,  respectively,  or  to  the 
free-space  field,  described  by  an  infinite,  denumerable  set  of 
modes.  We  will  refer  to  the  two-mode  coupling  as  Case  I,  and  to 
the  free-space  coupling  as  Case  ZZ,  distinguishing  between  these 
cases,  where  necessary,  by  the  subscripts  Z  and  ZZ,  respectively. 
In  the  absence  of  these  subscripts,  the  expression  is  applicable 
to  either  case.  For  simplicity,  we  consider  the  situation  in 
which  the  field  radiated  by  one  transition  has  a  negligible 
interaction  with  the  dipole  moment  of  the  other  transition;  in 
Case  I,  we  assume  that  the  “i3  mode  couples  only  to  the  <i,3> 
transition  and  in  Case  IZ,  we  assume  that  different  sub-sets  of 
modes  couple  to  the  dipole  moment  of  each  transition. 

Zn  boson  second-quantization  notation  til ,  the  atomic 
Hamiltonian  is  given  by  Hq  — Shu^ai"  oi,  and  the  dipole  moment  by 

-►  +4  +  + 

d  ■  Hi2«l«2  *  ^230203  *  h.c.,  where  ei  and  ai  are  the  usual 

annihilation  and  creation  operators  obeying  the  commutation  rules 
+ 

[ a i .  03J  ■  oij,  oil  other  commutators  vanishing.  Introducing  the 
reduced  variables  A3,  such  that  a]  ■  A3  and  using  the 

rotating  wave  approximation,  we  express  the  pump-atom  coupling  by 


K  I.jjAjA}  +  o23A2A*.1(“23-“>')  +  h.c.. 


where  ai3  *  -  (Eo  •  3  >  /<R  •  The  field  of  any  mode  can  be  described 

by  the  photon  annihilation  and  creation  operators  b  and  b+ ,  with 

*4* 

Cb(t>,  b  (t>3  *  1.  In  Case  I,  we  use  the  notation  b3j«,  with  the 

corresponding  reduced  variables  given  by  b3k  *  B3k  e“ia,3kt.  In 
Case  II,  we  use  the  notation  bj<,  with 
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b)<  ■  Bk«"iv  kt,  wher*  v  y,  is  the  frequency  of  the  k '  th  (free~space) 
mode.  Using,  again,  the  rotating  wave  approximation  £23  ,  we 
express  the  coupling  between  the  atom  and  the  field  by 
Haf  ■'h£AiA2$12  *  A2A^(2?  23^  *  h.c.. 


where 


riJ  BiJ(t) 


^i}  being  the  constant  that  describes  the  coupling  between  the 
<i,j)  transition  and  the  <i,j>  cavity  mode,  and 

Q±niM  *  J 

Yi3k  being  the  constant  that  describes  the  coupling  between  the 
<i,j)  transition  and  the  k'th  free-space  mode.  Introducing  the 
notation  i  =  (l/2)u>13-w  A  =<l/2>  <U12  “  w23>*  one  obtains,  as  atomic 
equations  of  motion 


A2  -  -l(o12A1.1<S+4)t+  0 23*3  e 


.  -i(5-A)t,  .  /D  JQ+  A  ) 

23A3  e  +  Ar-J12  +^23  A3  '  * 


A3  ■  -  i(o23A2e  +a2*023^  • 


Expression  for  the  field  in  terms  of  atomic  variables  are 
obtained  from  previously  derived  results.  These  are  133 

b«  <o>-1  iviji  YdtiAt<ti>AJ  ('I’'"'11  11  • i<3 
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where  ^j.3  la  the  loaa  conatant  of  the  (i,3>  mode  (energy  will 


decay  freely  as  exp  (-25j]t)),  and  [41 


-  lkij  A^(t)Aj(t),  i<j 


where  kj,3  ia  a  positive  constant  proportional  to  the  average 
value  (averaged  over  the  index  k>  of  lYufel2  in  the  neighborhood 
of  vj<  *  ui3#  and^^ii  ia  the  unperturbed  ,  that  is,  the 
operator  obtained  by  replacing  B)<(t)  by  B)<(0)  in  the  expression 


for  »J  i  jii  ( t)  .  (It  is  assumed  that  the  pump  is  turned  on  at  t  *  0.) 
We  will  consider,  presently,  only  the  situation  in  which  the 

loss  in  the  cavity-modes  is  sufficiently  large  so  that  the  field 
follows  the  atomic  polarization  adiabatically ;  in  other  words, 
j  is  larger  than  the  rate  of  change  of  In  this 

situation  one  can  write,  for  t>> , 


-0iSi-iki3IAt(t)V°  ' 


where 


kiJI  “  Kjl  , 


ijl  Tij  “ij 


Since  Q  l>  ■  0  for  both  Case  I  and  Case  II,  the  expression 
for  ^13 (t)  has  the  same  form  in  both  cases.  It  is  therefore 
unnecessary  to  treat  them  separately,  and  the  subscripts  I  and  II 
will  be  dropped  with  the  understanding  that  the  discussion. 


henceforth  refers  to  both  cases. 


Our  present  interest  lies  in  the  steady - state  atomic  level 
populations  and  dipole  moments.  We  need,  therefore,  equations 
for  the  expectation  values  of  all  bilinear  products  AiA3  with 
respect  to  an  (initially)  unexcited  radiation  field.  (It  shou 1 d 
be  noted  that  the  quantities  <A^Ai>  are  the  elements  Pj.3  of  the 
atomic  density  matrix  in  the  interaction  representation.)  The 
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derivation  of  such  equations  from  thoae  for  the  Ai'a  in  a 


convenient  form  is  outlined  briefly  in  the  following. 

We  first  write  expressions  for  all  derivatives  (d/dt)<A^A}> 
and  order  each  term  so  that  appears  on  the  extreme  left  and 

appears  on  the  extreme  right.  For  the  case  of  an  initially 
unexcited  radiation  field,  the^i^  term  refers  to  the  vacuum 
field,  and  the  kij  term  accounts  for  the  radiation  reaction. 
Substituting  for  each$i3<t>  these  two  terms,  and  utilizing  the 
relationship  < I  ««>  <t>  l>  ■  0,  we  drop  all 

explicit  vacuum-field  references.  The  above  ordering  of  the 
U'b  also  determines  the  ordering  of  the  four  atomic  factors  in 
the  ki3  terms,  which  we  now  rewrite  in  normal  order,  using  the 
appropriate  commutation  relationships.  Utilizing  another 
previously  derived  result  til  ,  namely, 

Aj,Ajl>  «  O  when  the  state  l>  describes  a  single  atom,  we  drop 
all  four-factor  terms  (remaining  after  the  normal  ordering). 
Finally,  the  substitutions 

Ajl  -  x^e-16*  ,  A 2  -  *2eiAt  *  A3  *  *3el6t  * 


Pij  m  <xixl  +  xtxj >  *  SiJ  "  i<xiVXiV 


ore  used  to  obtain  equations  of  motion  with  constant  coefficients 
for  nine  real  variables: 
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“l  “  °12S12  +  2k12n2  * 

“2  *  '  °12S12  +  °23S23  "  2k12n2  +  2k23n3  * 

n3  "  “  a23S23  "  2k23n3  * 

Sl2  -  -  («+A)P12  -  2c12(n1-n2)  -  °23P13-k12Si2 
S23  “  -  U-A)P23  -  2a23(n2-n3)  +  «12P13  “  (ki2+k23)S23  f 
§13  -  -  25P13  -  a23P12  +  a12P23  "  k23S13  * 

*12  "  (6+A)S12  +  °23S13  "  k12?12  , 

*23  "  (6-a)S23“  a12S13"  ^k12  +  k23)P23  ’ 

*13  “  26S13  +  °23S12  "  °12S23  “  k23P13  * 

These  equations  exhibit  2k ^ j  and  2a as  the  decay  rate  and  Rabl 
frequency,  respectively,  associated  with  the  <i,3>  transition. 
When  multiplied  by  the  appropriate  matrix  elements,  P12  and  P23 
describe  the  components  of  the  dipole-moment  expectation  values 
in  phase  with  the  field,  while  S12  and  S23  describe  those  in 
quadrature  with  the  field  (which  are  responsible  for  the  power 
absorption) .  It  should  be  noted  that  these  expectation  values 
represent  only  the  coherent  part  of  the  dipole  moment.  The 
incoherent  part  must  be  obtained  from  other  than  the  lowest 
moments.  The  coherence  of  the  radiation  will  be  discussed  later. 

III.  SOLUTION  OF  EQUATIONS  OF  MOTION 

As  mentioned  previously,  our  present  interest  lies  in  a 
steady-state  solution.  Setting  the  time  derivatives  equal  to 
zero,  and  replacing  one  of  the  first  three  equations,  (say,  that 
for  n2>  by  Zni  ■  1,  we  obtain  a  set  of  nine  inhomogeneous 
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equations  for  nine  unknowns.  An  exact  analytic  solution  la 
formally  possible  but  complicated,  and  we  resort  at  first  to  a 
numerical  solution.  For  this  purpose  all  the  constants,  which 
ore  inverses  of  the  time,  can  be  normalized  by  taking  a  12  to  be 
unity;  all  frequencies  and  decay  constants  are  then  given  in 
terms  of  the  <1,2>  Rabi  frequency.  Figure  1  shows  a 
three-dimensional  graph  of  n2/ni  vs.  fi  and  A.  The  resonance 
that  appears  in  n2/nj  for  large  A  is  striking,  and  motivates  an 
approximation  that  makes  analytic  solutions  of  the  steady-state 
equations  much  simpler.  We  consider  the  case  A>>a^3>>ki^.  If 
the  time  scale  is  chosen  so  that  ai3-l,  one  can  regard  A and 
kj.3  as  small  quantities  (of  first  order)  compared  to  unity,  and 
examine  the  solution  for  small  6.  The  four  equations  containing  A 
can  now  be  approximated  by  retaining  only  the  lowest  order  terms, 
that  is.  by  dropping  the  terms  (in  these  equations  only) 
containing  6  and  k.  The  resulting  set  of  steady  state  equations 
yields  a  solution  which  can  be  written  most  succintly  as  follows: 
let 


r  =  k 


12 


-  Q^2)/2A  , 


and 


R 


~r~r  'k23  +  “{-V2i  ’ 

°12  °23 
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With  this  notation,  we  have 
n2  -  [l+r(2+R)]"1  , 

n1  *  r(l+R)n2  , 


Figure  2  displays  both  the  numerical  solution  of  the  exact 
equations  and  the  analytic  solution  of  the  approximate  equations 
for  the  level-populations  as  functions  of  * .  It  is  seen  that, 
for  the  parameters  used,  the  approximation  is  good. 

IV.  DISCUSSION  OF  SOLUTIONS 

Since  the  quantity  F  reaches  a  minimum  at  <5  ■  6q»  the 
analytic  solutions  for  n2  and  n3,  as  well  as  those  for 
ISi2'>  1 523 1  ®nd  I  S 1 3 1  »  all  exhibit  resonant-type  maxima  as  a 
function  of  pump  frequency  at  <l/2)  ui3*  ^o*  that  is,  when  the 
pump  is  detuned  from  exact  two-photon  resonance  by  $o»  The 
half-maximum  resonance  width  about  for  either  n2  or  n3  is 


given  by 


lO 


and  that  for  n2'ni  ia  glvan  by 


2<6-S0>. 

u  JSoa» 


a12°23  .  .  2 
— r  +k23 

A 


Both  widths  ora  an  order  of  sag  nituda  smaller  than  the  Robi 

frequencies  2<*i3*  Population  inversion  is  achieved  either  in  the 

<2, 3)  transition  for  r>l,  or  in  the  <1,2>  transition  for 

2  2  2  2 

r<[l  *  <A  K23/®12  “23>J_i«  While  the  <2,3>  inversion  ratio  is 
frequency  independent  in  the  solution  of  the  approximate 
equations*  the  <1»2>  inversion  ratio  is  sharply  resonant*  and  is 
given  at  resonance  by 

(n,/n,)  -  r”1  (  1+A2k  2  a  “2  al?  )  ~1 

1‘  1  max  \  23  12  23/ 


The  transition  at  which  inversion  occurs  will  have  either  a 
higher  or  a  lower  frequency  than  the  pump  frequency*  the 
difference  being  approximately  A  •  It  is  also  interesting  to 
note  that  P3.3  goes  through  zero  at  resonance. 

It  is  clear  that  the  present  effect  may  have  interesting 
applications,  firstly,  because  population  inversion  for  a  given 
pair  of  levels  may  not  be  otherwise  achievable*  secondly*  because 
inversion  can  be  echieved  at  a  higher  frequency  than  the  pump 
frequency,  and  thirdly,  because  a  sharp  resonance  has  many 
obvious  uses  in  physics.  Discussion  of  the  details  of  such 
applications  is  beyond  the  scope  of  the  present  paper. 
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V.  RATE-PROCESSES  ANALYSIS 


It  is  desirable  to  offer  a  sore  intuitive  type  of 
explanation  of  the  present  resonance  effect  then  that 
contained  in  the  mathematics.  The  first  part  of  such  an 
explanation  is  the  observation  that  (l/2>ui3-$o  i»  Dust  the  pusp 
frequency  at  which  a  three-level  system  such  as  the  present*  but 
without  radiation  losses  (that  is,  without  coupling  to  the 
radiation  field),  behoves  like  a  two- level  system  for  A  >><*iDC5). 

Its  population,  if  intially  in  the  ground  state,  mainly 
oscillates  between  the  first  and  the  third  level  (a  Rabi-type 
oscillation)  with  frequency  2ai2a23/A*  while  the  second-level 
population  remains  small  and  oscillates  with  relatively 

high  f requency <-A ) E6D .  We  can  regard  our  present  radiative  system 
responding  to  the  pump  like  this  lossless  system,  pumped  from 
level  1  to  level  3  by  a  two-photon  process.  The  second  part  of 
the  explanation  is  the  observation  that  the  present  system 
behaves  like  a  three-level  system  with  respect  to  radiation 
losses,  relaxing  at  both  intermediate  frequencies  with  a 
consequent  distribution  of  population  to  all  levels. 

Since  the  above  explanation  invokes,  implicity,  rate 
processes,  it  is  of  interest  to  see  how  it  can  be  presented  in 
quantitative  form  using  rate  equations.  (It  should  be  noted  that 
the  explanation  is  applicable  only  to  the  range  of  parameter 
values  for  which  the  approximation  was  made.)  A  set  of  rate 
equations  for  three  levels  reads. 
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r»i  *  Pl2<n2~ni>  ♦  Ri3<n3-ni>  ♦  2ki2n2» 

r»2  *  Rl2<ni~n2>  *  R23(n3“n2>  *  2k23n3-2ki2n2 » 

113  *  R23<n2-n3>  *  Rl3<ni-n3>  -  2k23n3, 

wh«r«  Ri3  is  the  pumping  rate  between  levels  i  and  3.  In  order 
to  obtain  the  pumping  rate  between  consecutive  levels,  we  look  at 
the  equations  of  motion  for  a  two-level  system.  These  can  be 
obtained  immediately  from  the  three-level  equations  of  motion  by 
eliminating  all  references  to  either  the  first  or  the  third 
level.  Dropping  all  references  to  the  third  level,  we  obtain 

ni  *  a125i2  *  2ki2n2.  n 2  ■  “  ni, 

S12  •  -  <  6*  A)Pi2 -2a12(n1-n2)  -  k^S^, 

Pl2  ■  <  6  *  A  >Si2  “  *12Pl2» 

where  the  deviation  of  the  pump  frequency  from  resonance  is  now 
given  by  6 ♦  A.  Since  the  rate  equations  for  two  levels  reduce  to 

ni  »  Rl2<n2~nl>  *  2ki2n2#  n2  *  -  ril, 

we  must  have 


Rl2<r»2-ni>  *  “12512. 

Implicit  in  the  use  of  rate  equations  is  the  assumption  that 
rates  are  independent  of  the  state  of  the  system.  We  consider 
the  steady  state,  set  3i2  *  Pl2  *  0,  eliminate  P12  from  the  last 


two  equations  of  motion,  and  obtain 


Q12S12 


2a12kl2 
k^2  +  (6+A) 


2  (^2”“  1^  * 


We  return,  now,  to  the  three-level  rete  equation  and  compare 
them  with  the  three-level  equationa  of  motion,  obtaining  thua 

<212512  *  Rl2<"2“nl>  *  Rl3<n3-ni>. 

For  the  range  of  parametera  that  led  to  the  approximate  equations 
of  motion,  R12  is  ®  small  quantity  of  the  order  of  ki2/^2#  that 
is,  of  third  order.  (It  should  be  recalled  that  ku  and  A"* 
are  small  quantities  of  first  order,  and  (*13  is  of  order  unity.) 

From  the  steady  state  equations,  we  see  that  °^2Si2  must 
be  of  order  kj.2»  since  it  equals  **2ki2n2*  and  ni  is,  of  course, 
of  order  unity.  R12  is  therefore  negligible  compared  to  Ri3,  and 
R 1 3  must  be  of  the  same  order  as  Si2«  A  similar  argument  applies 
to  R23  relative  to  R13.  One  can  say  that  the  approximations  made 
in  the  equations  of  motion  find  their  counterpart,  in  the  rate 
equations,  in  the  neglect  of  R12  and  R23  compared  to  R13. 

Writing  the  three-level  rate  equations  in  accordance  with  this 
prescription  we  have 

hi  ?  Rl3<n3“ni>  *  2ki2”2* 

«2  =  2k23"3  “  2ki2n2» 

”3  -  Rl3<ni~n3)  ~  2k23H3. 

This  is  just  the  quantitative  form  of  the  qualitative  explanation 
offered  earlier.  If  we  now  set 


2°12°23k23 

'13  “  a2  .  2  .  ...  .  ,2 

A  k23  +  A (6-60) 


* 


these  rate  equations  will  yield  exactly  the  same  expression  for 
the  steady-state  populations  as  the  analytic  solution  of  the  full 
(approximate)  equations  of  motion.  It  should  be  noted  that. 
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while  rates  for  rate-aquation  purposes  are  usually  taken  to  be 
the  transition  probabilities  derived  by  perturbation  theory*  the 
above  expression  for  Rj.3*  with  the  detuning*  can  come  only  from 
consideration  of  the  full  set  of  coherent  equations.  Lastly*  as 
a  check*  we  compare  the  rate  of  energy  absorption  from  the  pump 
with  the  rate  of  energy  loss  to  the  field.  The  former  is  given 
by 

/h(«ui2  *  (*)23>Ri3<ni-n3> 
while  the  latter  is  given  by 
2/h  Wi2hl2n2  ♦  2ftw23*<23n3 

Substitution  of  the  above  expression  for  R13  and  the  analytic 
solutions  for  ni  in  these  two  rates  exhibits  their  equality. 

VI.  COHERENCE 

It  is  of  interest  to  examine  the  coherence  properties  of  the 
radiation  emitted  by  the  pumped  atom.  For  simplicity*  we  use  the 
notation  applicable  to  Case  I*  in  which  the  atom  radiates  into 
two  lossy  cavity  modes.  The  (complex)  field  operator  for  the 
<l;j>  mode  la  proportional  to  Bu  expC-i^ut).  By  factoring  out 
"*13  from  the  expression  for  ill(t)  in  Sec.  II*  we  obtain 

ivr  + 

B^U)  -  (0)  -  Ai(t)Aj(t)  ,  i  <  j  , 

where*  it  is  recalled.  *  0.  The  coherent  field  is 

proportional  to  <Bj.j  (t)  >exp<-iu>i3t)  .  The  expectation  value  of 
the  total  energy  in  the  <i,j)  mode  is  'fiwu  OijBij),  while  the 
energy  associated  with  the  coherent  field  is  HSu)  i3<B*i3>  <Bi3>. 
From  the  above  operator  expression  for  Bij,  we  have 


The  order  of  Magnitude  of  this  quantity  is  the  same  as  that  of 
<kij/<»iO)2  *  <<*i3/A>2,  which  is  a  small  quantity  of  second  order 
for  the  parameter  range  under  consideration.  We  see.  therefore, 
that  most  of  the  radiation  is  incoherent.  The  frequency  of  the 
coherent  radiation  is  that  of  the  pump,  as  is  to  be  expected;  the 
incoherent  radiation,  on  the  other  hand,  may  be  expected  to  lie, 
largely,  in  the  frequency  regions  near  the  two  transition 
frequencies.  A  quantitative  discussion  of  the  spectral 
distribution  of  the  radiation  is  beyond  the  scope  of  the  present 
paper . 

VII.  CONCLUSION  AND  SUMMARY 

Single-frequency  two-photon  pumping  of  a  three-level  system 
that  radiates  naturally  only  at  the  two  intermediate  frequencies 
has  been  considered.  It  has  been  shown  that  when  the  pump  is 
slightly  detuned  from  two-photon  resonance  and  not  too  close  to 
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one-photon  resonance  for  one  of  the  Intermediate  frequencies, 
population  inversion  con  be  obtsined  either  in  the  upper  or  in 
the  lower  psir  of  levels,  depending  mainly  on  the  ratio  of  the 
two  radiative  decay  constants.  This  effect  furnishes  a  novel 
method  of  population  inversion  for  laser  purposes  and  a  offers 
the  possibility  of  obtaining  population  inversion  in  a  transition 
with  a  frequency  higher  than  the  pump  frequency. 
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Fig.  1 


The  steady  state  population  ratio  in  a  three-level 

system  as  a  function  of  and  A  =  H (‘‘,i2_w23^  for 


1,  o. 


2,  k12  ■  0.01  and  k23  »  0.04. 


choose  (*23^12^  =  ^°23^a12^  ) 
aij-| j |  and  kAj  «.|wij|2.) 


since 
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A  comparison  of  the  exact  numerical  solution  (solid 
lines)  and  approximate  analytic  expressions  (dashed 
lines)  for  the  steady  state  populations  as  functions 
of  6,  with  a.-  =  1»  a-~  = 


2,  k. -  =  0.01,  k?_ 


=  0.04  and 
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